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نشربه بوم شناسی کشاورزی 


9 ۷ :۳۱08260226 
که 2 ۲ ۳ 


مقایسه روش‌های زمین آمار و سنحش از دور در برآورد عملکرد دانه برنج (1,۰ 3۹0/۳0۵ 0#)۳۳۳2۵) 
در استان گیلان 


پویا اعلایی باز کیایی . بهنام کامکار"" ابراهیم امیری ؛ حسین کاظمی * و مجتبی رضایی" 
تاریخ دریافت: ۱۳۹۹/۰۶/۰۷ 
تاریخ بازنگری: ۱۴۰۰/۰۳/۲۷ 


تاریخ پذیرش: ۱۴۰۰/۰۳/۳۱ 


اعلایی بازکیایی» پ.» کامکاره با.» آمیری» ۷ کاظمی. .۰ و ر ضایی. م۰ ۳-۳۰۱ مقاد سه روش‌های زمین‌آمار ۴ سنحش از دور در برآورد عملکرد 
دانه برنج (م1 ۶0۲۵ 0:2۵) در استان گیلان. بوم‌شناسی کشاورزی ۳(۱۴): ۰۵۷۹-۵۹۹ 


چکیده 


به‌منظور مقایسه روش‌های مبتنی بر سنجش از دور و زمین‌آمار برای برآورد عملکرد برنج (.1 90/۷۵ 07(20) اراضی کشت برنج استان گیلان, 
پژوهشی در سال‌های زراعی ۱۳۹۵ و ۱۳۹۶ انجام شد. جهت انجام عملیات میدانی» ۳۲۰ نمونه عملکرد دانه در مرحله رسیدگی فیزیولوژیکی گیاه برنج از 
۸ زار هکتار اراضی برنج برداشت شد. در این تحقیق. از ۳۳ روش-مدل زمین‌امار برای درون‌یابی مقدار عملکرد دانه استفاده شد و سپس دقت 
روش‌های درون‌یابی با کمک معیارهای آماری مختلف ارزیابی شد. در روش مبتنی بر سنجش از دور از تصاویر سنجنده تصویربردار عملیاتی زمین 
تما ماهواره اقدمت‌سارو تضاویر تیتجنده کصویزیرخار ماهوازه متتیل۲ امتفاده شد.‌همت قراخ گیاهی ام تصاویر انتخراج و ارباط بیخ آن‌ها چ 
متغیر عملکرد با | ستفاده از رابطه رگر سیونی | ستخراج و نق شه عملکرد با کمک صاویر تهیه و ارزیابی شد. نتایج ارزیابی روش‌های درون‌یابی بیان گر 
برتری روش-مدل کریجینگ معمولی - ٩18016‏ نسبت به دیگر مدل‌ها بود. در اين مطالعه. شاخص ]1۷ در بین شاخص‌های گیاهی بیش‌ترین دقت را 
در پیش‌بینی عملکرد ارا ضی برنج دا شته است. مقایسه روش زمین‌آمار و روش مبتنی بر سنجش از دور بر اساس ضریب تبیین و جذر میانگین مربعات 
خطا در برآورد عملکرد دانه در سطح | ستان گیلان ذشان داد که هر دو روش دقت قابل قبولی دا شتند. امّا با توجه به توانایی سنجش از دور در تفکیک 
نقطه‌ای بازتاب نوری پدیده‌ها و پیش‌بینی عملکرد با تفکیک مکانی بالاه در این روش دقت بالاتری در تخمین عملکرد حاصل شد. 


واژه‌های کلیدی: سنتینل- ۲ عملکرد واقعی» کریجینگ» لندست- ۸ 


۱- دکتری زراعت. دانشگاه علوم کشاورزی و منابع طبیعی گرگان. ایران. 
۲- استاد. گروه اگروتکنولوژی. دانشکده کشاورزی» دانشگاه فردوسی مشهد. مشهد. ایران. 
۲- استاد» گروه مهندسی آب. دانشکده فنی. دانشگاه آزاد اسلامیء واحد لاهیجان؛ ایران. 
۴- دانشیار. گروه زراعت. دانشگاه علوم کشاورزی و منابع طبیعی گرگان» ایران. 
۵- استادیا, موسسه تحقیقات برنج کشور سازمان تحقیقات. آموزش و ترویج کشاورزی» رشت. ایران. 
(#- نویسنده مسئول: 1600۷۹۱۱ ۵ ۲ هصق :اتفصط) 
308.0- 1022067/22 :10۵۲ 


2۸۷۲ نشر به بوم شناسی کشاورزی. جلد و۹ شماره۲ . باییز ۱۴۰۱ 


مقدمه 

برنج (. 501/۷۵ 012۵) دومین محصول زراعی از نظر سطح زیر 
کشت و قاسینکننده. غنای, بیش از تیضی از جفعیت جهان استة 
(2019 ,2 61 21۳21؟1). یک‌پنحم جمعیت جهان يا بیش از یک 
میلیارد خانواده در آسیاء آفریقا و آمریکا به‌عنوان منابع اصلی اشتغال و 
معیشت خود به نظام‌های تولید برنج وابسته هستند. همچنین برنج در 
مبارزه با گرسنگی و فقر جهانی در خط مقدم قرار دارد 6 /۳۵۳۵) 
سال ۲۰۱۸ دارای ۵۸۰ هزار هکتار سطح زیر کشت و ۱/۹۹ میلیون تن 
تولید برنج بوده است. واردات برنج در ايران در سال ۲۰۱۸ حدود ۱/۲ 
میلیون تن بوده است (2018 ,۳۴۸). این نتایج دلالت بر آن دارد که 
تولید داخلی کفایت تأمین نیاز کشور را نداشته و بخشی از نیاز برنج 
مصرفی از طریق واردات تأمین شده است. 

با توجه به مطالب ذکر شده درباره گیاه برنج نظارت بر تولید آن در 
سطح وسیع امری مهم است (2019 ,.21 61 ۷۲2). مقایسه روش‌های 
مختلف بررسی عملکرد گیاهان زراعی در سطح وسیح نشان داد که 
روش‌های سنتی با انتخاب تصادفی اراضیء با توجه به وجود تغییرات 
جغرافیایی» اقلیمی و خاکیء اطلاعات مناسبی در مورد میزان عملکرد 
دهد (2013 بللعاما :2019 له 6 ۷۲2۵). در تعیین عملکرد 
محصولات مختلف» سنجش از دور نقش مهمی در ارزیابی مکانی و 
زمانی عوامل اقلیمی» خاک» عوامل زنده و روش‌های مدیریتی اثرگذار 
بر عملکرد دارد که مهم‌ترین مزیت آن پردازش داده‌های مکانی وسیح 
می‌باشد (2020 ,.له 6 1601060۲01 طع۵11220). سنجش از دور را 
می‌توان علم کسب اطلاعات از سطح زمین بدون تماس فیزیکی با آن 
امید که از جمله استفاده از عکس‌های هوایی 9 تصاویر ماهواره‌ای 


می‌باشد. روش سنجش از دور عبارت است از اندازه‌گیری و ثبت آنرژی 


بازتابی یا متشر شده الکترومفناطیسی از سح زمین و جو که در بین 
پدیده‌های مختلف دارای میزان متفاوتی است (2014 ,2057ظ). مزیت 
عمده تصاویر ماهواره‌ای قابلیت استفاده از آن به‌صورت رقومی است که 
می‌تواند به‌وسیله رایانه‌ها تجزیه و تحلیل شده و در زمان‌های مختلف 
در دسترس قرار گیرند (2006 ,.21 66 0002771)). یک راه مناسب 
بهره‌گیری از داده‌های تصاویر ماهواره‌ای استخراج شاخص‌های گیاهی 
مناسب با منطقه و پوشش زمین است (2014 ,20527ظ). شاخص‌های 
گیاهی تبدیل‌های ریاضی هستند که بر اساس باندهای مختلف سنجنده 
تعریف شده و برای ارزیابی گیاهان در تصاویر چندطیفی ماهواره‌ای 
طراحی شده‌اند (2014 ,320527). از جمله شاخص‌های گیاهی 
پرکاربرد می‌توان به ۱۲۲۷۲۳ و "1۷۲ اشاره نمود (2014 ,۱20527). 
زار کو-تجادا و همکاران (2005 ببلة 66 22۲00-16[202) مشاهده 
کردند که شاخص‌های "لها 19۷۲ ۱۲۳۷۴ و 05۸۷۲ با 
استفاده از ماهواره 110705 با وضوح مکانی چهار متر بهترین ارتباط 
را با عملکرد پنبه (م] ۷۸ ۶۱رععم)) در کالیفرنیا داشته‌اند. 
مو و همکاران (2005 ,.21 66 ۷]0) با استفاده از مدل‌سازی به پیش‌بینی 
مناسپی از عملکرد گیاهان زراعی چهارکرینی و سه‌کرینی» کارایی 
مصرف آب و کارایی تبخیر و تعرق با استفاده از ۲۷1 حاصل از 
سنجنده 00۵۵-۸۷11817 رسیدند. دنت و همکاران .21 61 ۲۵۵1) 
(2008 گزارش کردند که ایجاد نقشه مناسب شاخص سطح برگ با 
تصاویر ماهواره‌ای به افزایش دقت برآورد عملکرد با استفاده از مدل 
0۳۳5-۷۷۵۱۳ کمک می‌کند. (2011 ماعع26:0 ک نط) از 
طریق برقراری ارتباط بین ۱۷۲۷۳ (حاصل از تصویر ماهواره -612) 
5 با گستره دید ۲۳۲۰ کیلومتری) و حداکثر ظرفیت کاتالیزوری 
روبیسکو ۷000 الگوی ظرفیت فتوسنتزی منطقه مورد مطالعه خود 
را به‌دست آوردند. محمدی احمدمحمودی و همکاران ن4صصهطم۷) 
(2015 ,.له 6 لامصطد]/ ۸۳۳۱20 مشاهده کردند که شاخص 
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۲ تصویر سنجنده "+7۲1 مربوط به اواخر مرحله پنجه‌زنی گندم 
(س ده 71:6۱ عملکرد مزرعه را با دقت پیش‌بینی می‌کند. 

زمین‌آمار یکی دیگر روش‌های بررسی تغبیرات عملکرد در سطوح 
وسیع بر اساس الگوریتم عددی صفات پیوسته است. کریجینگ یکی از 
روش‌های زمین‌آمار برای درون‌یابی داده‌ها بر اساس واریانس مکانی 
است (2015 ,.21 ]6 01امصطع]۱۷ 0عحصطم محصجطهط۷0). سامانه 
اطلاعات جغرافیایی در برنامه‌ریزی مکانی با تعیین پراکنش پدیده‌ها و 
روی هم‌گذاری نقشه‌ها و تفسیر داده‌های بوم‌شناختی» در مراحل 
مختلف برنامه‌ریزی کاربرد دارد (2005 ,5ع61 > 13251606). بررسی 
تنوع زیستی علف‌های هرز در سامانه‌های زراعی تحت کشت گندم 
شهرستان گرگان (2011 ,010256۳01))» بررسی رابطه شوری با عملکرد 
گندم در مزرعه نمونه ارتش گرگان (2019 ,.21 6 ت۸2 بررسی 
رابطه عناصر غذایی خاک با عملکرد گندم ,۸20 طلع5ه۳0 [2۵6ع۳5) 
(2013 برآورد خلاء عملکرد در حوضه گرگان با سامانه اطلاعات 
جغرافیایی و مدل ۹5۷۲ (2014 ,20527ظ) و پایش تولید زیست‌توده 
در مزارع گندم ارتش استان کلستان 24صطه ن4مصصحطه]) 
(2015 ,یله ]6 ۱۷2۳00101 از جمله موارد کاربرد سامانه اطلاعات 
جغرافیایی در علوم کشاورزی می‌باشند. 

بنا بر مطالعات ذکر شده سنجش از دور و روش زمین‌آمار رهیافت- 
های کاربردی در زمینه بررسی عملکرد و متغیرهای گیاهی در سطح 
وسیع می‌باشند. با توجه به اهمیتی که تولید برنج در تأمین امنیت غذایی 
برنج کشور درد بررسی میزان عملکرد بهدست آمده در سطح استان 
گیلان که دومین استان از نظر سطح زیر کشت برنج (۲۳۸ هزار هکتار) 
در کشور می‌باشد. می‌تواند اطلاعات مناسبی در پیشبرد اهدافی چون 
بررسی میزان خلاء عملکرد در این منطقه فراهم نماید؛ بنابراین اين 
پژوهش با هدف بررسی عملکرد برنج با استفاده از دو روش زمین‌آمار و 
رهیافت سنجش از دور در سال‌های ۱۳۹۵ و ۱۳۹۶ در سطح استان 
گیلان انجام شد. 


مواد و روش‌ها 
منطقه مورد مطالعه 

این تحقیق در سال‌های زراعی ۱۳۹۵ و ۱۳۹۶ در اراضی شالیزاری 
استان گیلان انجام شد. استان گیلان دارای ۱۶ شهرستان با ارتفاع 
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۱- تا ۲۶۸۷ متر از سطح دریا می‌باشد که دشت‌های آبرفتی و مناطق 
کوهستانی از سمت شمال به جنوب امتداد یافته است (سازمان مدیریت 
و برنامه‌ریزی استان گیلان). استان گیلان در مختصات جغرافیایی ۲۶ 
درجه و ۲۷ دقیقه تا ۲۸ درجه و ۲۷ دقیقه عرض شمالی و ۴۸ درجه و 
۲ دقيقه تا ۴۹ درجه و ۴۵ دقیقه طول شرقی قرار دارد. این استان» بر 
اساس طبقه‌بندی کوین-گیگر (بر اساس داده ۲۵ ساله) دارای اقلیم 
تا برداشت برنج در این استان از اوایل فروردین (خانه‌گیری نش) تا پایان 
شهریور (برداشت محصول) می‌باشد. وسعت اراضی شالیزاری استان 
علی کاظمی» جمشیدجو» خزر بیش‌ترین سطح زیر کشت را داشتند 
(2016 مصهاتیدن طا ممتاهم2نصهع:0 120 -عتیمآنه‌توش). 


تعریف نقاط نمونه‌برداری 

منطقه به‌وسیله دستگاه 0۳5 مدل ۵۵۰ ص02 ثبت شد. سپس 
کمینه تعداد شالیزارهای مورد نیاز برای بررسی به‌صورت تصادفی و بر 
اساس روش تعیین تعداد نمونه کوکران در سطح استان تعیین و جامعه 
آماری (۳۲۰ مزرعه) انتخاب شد (2019 ,.ه !6 صمن‌مطستم۳). در 
زمان برداشت محصول به‌منظور تعیین عملکرد دانه. به شالیزارهای 
پیش از ریزش دانه با کد زادوکس )٩۲-"‏ با همکاری ناظرین جهاد 
کشاورزی صورت گرفت. سپس نمونه‌های گیاهی به آزمایشگاه منتقل 
و پس از تعیین درصد رطوبت» میزان عملکرد دانه با رطوبت ۱۴ درصد 
محاسبه شد. مختصات جغرافیایی مزارع در شکل (۱) ارائه شده است. 
از بین نقاط نمونه‌برداری» ۷۰ درصد آن برای برآورد عملکرد دانه برنج 
با روش سنجش از دور و زمین‌آمار در سطح استان گیلان و ۳۰ درصد 
باقی‌مانده برای صحت‌سنجی روش‌های مذکور استفاده شد ,20527ظ) 
2014 


روش درون‌یابی 
در این پژوهش» روش درون‌یابی کریجینگ به‌دلیل دقت بالا و در 
نظر گرفتن ساختار مکانی و همبستگی بالا در داده‌ها مورد استفاده 
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قرار گرفت (2015 ,.له 6 ن4تامصطع]۷ معط تمصصه‌طهما) 
بدین منظور آزمون ران تست" جهت اطمینان از همگنی داده‌ها انجام 
نرمال بودن توزیع داده‌ها با استفاده از شاخص‌های مرکزی شامل میانه 
و میانگین و شاخص‌های پراکندگی شامل چولگی و کشیدگی مورد 
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شکل ۱- موقعیت جغرافیایی اراضی سالیزاری و مزارع نمونه‌برداری در استان گیلان در سال ۱۳۹۵ (الف) و ۱۳۹۲ (ب) 
() 2017 200 (2) 2016 1 0۲۵۷۱۵۸۵ حهآنبی) صا فص رو ۱60مرصعی 0صه صمتاه یت م۲۳۱6 ما مها ۵۶ ممتای‌م۱ آمعنطام‌جدوهعی - 1 ۲۱۵۰ 


مبنای نرمال بودن داده‌ها به‌دست آمدن مقدار چولگی و کشیدگی 
در محدوده ۳+ تا ۳- بود (2014 ,2052ظ). به‌منظور بررسی ساختار 
مکانی داده‌ها نیم‌تغییرنما مورد بررسی قرار گرفت که ساختار ارتباط 
مکانی بین نمونه‌ها را نشان می‌دهد (2001 ,013۷۵۲ > تعادداه ۷۷). 
متغیرهای نیم‌تغییرنما شامل اثر قطعه‌ای " (خطاهای اندازه‌گیری یا منابع 
تغییر مکانی در فواصل کوچک‌تر» دامنه یا عاع تأثیر" (ملفهای از 
مدل کوواریانس یا نیم‌تغبیرنما که فاصله‌ای را تشان می‌دهد که فراتر 
از ان شمیشگی ود کادییی عی‌ها کم است با اضلا پحوه بارداع 
آستانه يا سقف" (مقدار به‌دست آمده در مدل نیم‌تغییرنما در محور ۷ در 
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نقطه شعاع تأثیر) می‌باشند (2000 ,.21 6۱ 1[1561). برای بررسی ساختار 
مکانی داده‌های مورد استفاده نسبت اثر قطعه‌ای به سقف در هر نیم- 
تغییرنما بررسی شد. همچنین مقدار عددی شعاع تأثیر* بر مبنای متوسط 
فاصله نقاط از هم و با استفاده از تابع نزدیک‌ترین مجاور محاسبه و 
استفاده شد (2014 ,20527). تعداد گام و نوع شکل تخمین" نیز با 
آزمون و خطا بر اساس کم‌ترین ريشه دوم میانگین مربعات خطا 
([۹)2۷]5 با حد خطای ۲۰ درصد میانگین تعیین شد. سپس ایستگاه- 
های دارای خطای بالا مشخص و حذف شدند و بار دیگر میان‌یابی با 
تعیین شعاع تأثیره تعداد گام و مناسب‌ترین شکل تخمین انجام شد. 
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به‌منظور اطمینان از نرمال بودن داده‌هاء توزیع فراوانی آن‌ها با کمک 
چولگی و کشیدگی توسط نرم‌افزار آماری 5۸5 و رویه 2۲0۵ 
عوت2 ۷ مورد بررسی قرار گرفت (2006 ,226ع56 عک تصع/۳۲۵). 
پس از تأیید نرمال بودن داده‌هاء ۳۳ روش-مدل (سه روش کریجینگ 
معمولی» جهانی و ساده 9 ۱۱ مدل دایره‌ای 4 کروی * نمایی * گوسی * 
چهاروجهی کروی* پنج‌وجهی کروی* درجه دو منطقی" ثابت* اثر 
روزنه‌ای" 1۳65561 و 1-۳65501) مورد ارزیابی قرار گرفتند که برای 
هر متفیر میزان خطای موجود در تخمین نقاط مجهول مورد بررسی 
قرار گرفت. جهت تشخیص مناسب‌ترین روش؛ میزان 15 در هر 
روش محاسبه و ارزیابی شدند (2012 ,.21 6۲ 12276101). به‌طور معمول» 
هر قدر مقدار 5۷18۳ کمتر باشده صحت روش بیشتر است. از نظر 
تتوری هرگاه اين معیار برابر صفر شود نمایانگر این است که مقنار 
تخمین زده شده یک کمیت» دقیقاً برابر مقدار واقعی آن است 
(2012 .له )6 1>226۳01). 


تصاویر ماهواره‌ای 

در این پژوهش, از پنج تصویر چند طیفی سنجنده 01,1۲ ماهواره 
لندست-۸ که دارای نه باند بوده و به تاریخ‌های ۲۹ خرداده ۱٩‏ مرداد و 
۱ مرداد در سال ۱۳۹۵ و تاریخ‌های یک مرداد و ۱۷ مرداد در سال 
۶ تعلق داشتند. استفاده شد. تصاویر بر اساس درصد ابرناکی کمتر 
از دو درصد در محدوده زراعی برش داده شده استان انتخاب شدند. 
همچنین از تصاویر سنحنده ]1۷۲5 ماهواره سنتینل-۲ در تاریخ نهم تیر 
و ۲۲ شهریور در سال ۱۳۹۶ استفاده گردید. در این پژوهش. دقت 
حاصل از کاربرد تصاویر دو ماهواره لندست-۸ و سنتینل-۲ با قدرت 
تفکیک مکانی متفاوت در برآورد عملکرد دانه برنج مورد ارزیابی قرار 
گرفت. تصاویر دریافت شده نیاز به تصحیحات هندسی و اتمسفری 
نداشت. در این پژوهش» هشت شاخص گیاهی شامل ۰۲۷۲ ۷۲ 
۰۲9/۸۷۲ ۳۷ ]1۷۷ ]1:۷ و 1۷ مورد استفاده قرار 
گرفت (جدول ۱). باندهای مورد استفاده از هر سنجنده با ذکر محدوده 
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طول موج در جدول ۲ ارائه گردید. باندهای مورد استفاده از هر سنجنده 
شاخص‌های مذکور از نرم‌افزار 1201151 نسخه 18.08 ۳6۴۲ 
استفاده شد. سپس با کمک شاخص‌های گیاهی» نقشه عملکرد برنج 
تهیه و به‌منظور صحت‌سنجی نقشه عملکرد دانه برنج حاصل از 
شاخص‌های گیاهی, مقایسه عملکرد پیش‌بینی شده و عملکرد ثبت شده 
صورت گرفت. 

پس از محاسبه شاخص‌های گیاهی. مقادیر متناظر این شاخص‌ها 
در محل نقاط نمونهگیری استخراج و ارتباط بین این شاخص‌ها با 
عملکرد دانه برنج بررسی شد و معنی‌داری این شاخص‌ها با کمک 
رگرسیون ساده خطی مورد ارزیابی قرار گرفت و شاخص‌های مطلوب 
گیاهی با عملکرد دانه برآورد گردید و در محیط نرم‌افزار 10.3 ۸۵۲6615 
و با کمک رابطه ریاضی به‌دست آمده» عملکرد دانه از روی شاخص‌های 
گیاهی به‌دست آمده از تصاویر تخمین زده شدند. با استفاده از مقادیر 
به‌دست آمده عملکرد دانه به‌وسیله شاخص‌های مبتنی بر سنجش از 
دور و روش زمین‌آمان نقشه‌های عملکرد واقعی محدوده به‌دست آمد. 
سپس بهترین نقشه تخمین عملکرد با کمک ريشه دوم میانگین مربعات 
خطا و ضریب تبیبن (2) انتخاب 9 ارائه ۳ آزمون مورد استفاده برای 
بررسی معنی‌داری 22 آزمون ۳ در سطح احتمال پنج درصد بوده است. 
در این پژوهش از نرم‌افزارهای 10.3 ۵:۰015 5.1 سح 2161 
نسخه 18.08 1۳۳5۳۲ و نرم‌افزار آماری 9.4 5۸6 که6) 
(2015 ,6)ااتا105 استفاده گردید. 


بررسی عملکرد وأقعی با استفاده از فناوری سنجش از دور 

روابط خطی بین شاخص‌های گیاهی و عملکرد برنج به همراه 
ضریب تبیین و ضریب تغییرات برای سال‌های ۱۳۹۵ و ۱۳۹۶ در جدول 
(۴ و ۵) نشان داده شده است. در سال ۱۳۹۵ بر اساس مقادیر به‌دست 
آمده جذر میانگین مربعات خطا (جدول 4۴ شاخص ۷1 به‌دست آمده 


7- 2110821 60 

اطعا -5 

9- ۳۲۵18 1 

10- 00۵۵۲2010821 12۵0 ۳۵۵6۲ 


ور( نشر به بوم شناسی کشاورزی. جلد و۹ شماره۲ . باییز ۱۴۰۱ 


از تصاویر ماهواره لندست-۸ در تاریخ ۱۹ مرداد دارای کم‌ترین مقدار و 
ضریب تبیین معنی‌دار در روابط به‌دست آمده بود. در سال ۱۳۹۶ نیز بر 
اساس مقادیر به‌دست آمده جذر میانگین مربعات خطاء شاخص 1۷۲ 
به‌دست آمده از تصاویر ماهواره لندست-۸ در تاریخ ۱۷ مرداد کم‌ترین 
مقدار (۲۱۸:۰۳) و ضریب تبیین معنی‌دار را نشان داد. شاخص ۷ از 
نسبت نور مادون‌قرمز به نور قرمز بهدست می‌آید که در سطح کاملا 
صاف برابر صفر است که با رشد گیاه مقدار آن افزایش می‌یابد 1>10) 
(2009 ,211 ۷27 ۶. در این شاخص امواج الکترومغناطیسی منعکس 


شده تحت تأثیر ساختار هندسی تاج‌پوشش از جمله اندازه شکل و 
آرايش اجزای تاج‌پوشش است (2013 ,21 6 160206). در اين 
پژوهش, به نظر می‌رسد شاخص 18۷1 با توجه به خصوصیات ذکر شده 
در مورد این شاخص, با توصیف بهتر رشد گیاه برنج دقت بالاتری در 
برآورد عملکرد برنج در محدوده مورد مطالعه داشته است. در مطالعات 
دیگر از شاخص 18۷ جهت بررسی متفیرهای گیاهی» خاک و اقلیم 
استفاده شد. 


جدول ۱- روابط شاخص‌های گیاهی 
۰ ۲۸۱۵ 622)21100 ۷ -1 2016 1 


منبع رابطه علامت 
۱۹ ۳02 مار 
1۰ 6 01بامحطع۱۷ 0ححصطمه ت2صصصصهطام]1۷ وووطا - روط ۳۷1 شاخص گیاهی تبدیل شده 
(2015) 5 + (وموط + چرط) 6ص طمت۵ا6ع۷6 ۲285]01۲۲۵0 ]1 
۰ 6 زنامحصطه]۷ 0محصظم ته‌محصعصفطم] ۷( یط شاخص گیاهی نسبت 
1۷ 3 ۲ 
(2015) مصوط ۲ ۷6۵۵۵1108 20[ 
(2014) 1320527 برجم سس ] .ده ات وی 
باج+ومو۲ا روط 


۵9 61 21. )2018( 


60 61 21. )2018( 


(2015) 2821ظ » ق7016 


(ممعتط 
۰ 6 1بامحصطع۷( 0محصظه تلمحصصصعطام]۷[ مصود < روط 
(2015) مموط + ورد 
۰ 6 1بامحصطع]۱۷ 0محصظه تمصصحصع‌طم]1۷ 
(2015) حصوط تک ۲ 


(ظ - ووود < 2 - ورن <2 
2-0 مووطا ج ور < 2 
طا< وووت < 2 - ورد 


+ ۷1 
جورورط /(00 0 ورری() 


۲ ۷۵۵۵۵110۴ 0ع)ویا 5011-20 
شاخص گیاهی تغییریافته منظم شده خاک 


۷-+--1*( جمناهامعه۷ 0عامنازله]ز0و 0مصصممصهت 
1۴ 
شاخص گیاهی عمودی 
۳۷ ۴ 
۴ ۷۵۵۵۵0۵ ۳۲۵۱۵۵۵2 
شاخص تفاضل آبی نرمال شده 
۷( 
۷/۵۲۵۲ 01111016060 ۱۱۵۱۵۵11260( 

۷ شاخص تفاوت نرمال شده گیاهی 


۷۵۵۵۲20۵ 01111670060 2760ت1هص۱۱0 
ار« :۷۳۳ 
6 ۷6۵۵۱۵100 121116۲6066[ 


حصرورط: باند قرمز, 94 باند مادون‌قرمز نزدیک» جمعتن 0 باند سبزء ر1 عامل تنظیم برای کاهش اثر خاک روی پوشش گیاهی می‌باشد که دارای دو مقدار ۰/۵ و ۱ می‌باشد. ۰/۵ 


برای تراکم پوشش گیاهی متوسط و ۱ برای مناطق با پوشش گیاهی کم‌تراکم. عامل تصحیح خاک ۰/۵ در نظر گرفته شد (1988 ,عا306) 
طمنط محمتاهامعع۷ م۵ امد 0۶ )660 معط عمتا0ع۲ ما ما20 اصمصص منز :با یصوها همرت :مممین‌ها یصقها 40متتصة بقع تمیدها یصه‌ها مک تصصوعرط 
٩011 201(105)060 0۲‏ .راتعصعل حمتاهاعع۷6 10۷ 6 قدهتج ۵۶ 1 مه بواتعصعل ممتاهاععع۷ تلع ت۶0 0.5 :1 قصع 0.5 ۵۶ ومیاله۷ ۲۳۵ فقط 
19868(۰ ,۲۲۵666) 0.5 6088100760 ۲۷۷28 


جدول ۲- مشخصات تصاویر ماهواره‌ای دریافت شده 
٩261116 ۵۵65‏ ۲66۵۲۷۵0 01 ۵/215( -2 2016 1 


درصد ابرناکی ردیف مسیر تاریخ سنجنده ماهواره 
امه 0100۲0 10۲ ۱۳ ۱۱۳۹۹۵ 900۴ تناااتنن 
238 +034 166 ۹ خرداد (2016 ,18 06[) 01 لندست ۸ (8 0526ع) 
٩ 166 034 024‏ مرداد (2016 ,9 )5اعنظ) ]0 لندست ۸ (8 200526]) 
0 +034 166 ۱ مرداد (2016 ,21 )5اعظ) 99 لندست ۸ (8 1/200626) 
002 +034 166 ۱ مرداد (2017 ,23 17نا[) 017 لندست ۸ (8 0526صع) 
0 034 166 ۷ مرداد (2017 ,8 )5اعناض) ]01 لندست ۸ (8 100521) 
درصد ابرناکی شماره ردیف سطح تاریخ ماهواره 
010090 رباص ما1 ۳۱۵۵۲0 ۱۳3 عانلاها 
٩ 5۳۲1۳۲۳ 1399 0‏ تیر (2017 ,30 ۲086) ستینل-۲ (2 امصناجهو) 
٩۳۲] 39108 124‏ ۲۲ شهریور (2017 ,13 ]500 2-1 سنتینل-۲ (2 [5606) 


اعلایی و همکاران. مقایسه روش‌های زمین آمار و سنجش از دور در برآورد  ...‏ ۵۸۷ 


جدول ۳- مشخصات باندهای مورد استفاده در تصاویر ماهواره‌ای 
٩2161۱۱) ۷۲‏ صا 0و ول‌صهها ۵1 ممناوام)00۵2۲2 -3 م1201 


طول موج شماره باند 


طول موج شماره باند 


۰ ماهواره ۰ ماهواره باند 
۱۹3 ِِِ ِِ- (س منامام 1۵۲6۲۵8۵6۵6 ما ِِ 0 3( ۳90( 
سنتینل-۲ لندست ۸ ۱ 
0.533-0 3 9 0.533-0 3 5 
(2 6011061) (8 12020526) ات6 
متس لندست ۸ قرمز 
0( 0.636-3 4 ِ 0.636-3 4 گرعز 
(2 6011061) (8 1,290521) ۳20 
0851079 8 سنتینل ۲۳ 0851-09 5 ی 
(2 6011061) (8 2ع0صق) .۰ 1812760 


۳۱20000 _وتعو] ههد _قن]_5-1574اص1/واصمصنهم/وع1ی)/6ناهم/207/51)۵5 .1525 125://1200521 :2) 
01 1/567 _2-[1/5۵0060 1 9/247904/6852)صعصصیاع 16/0 عوم: امصتج1005://5۵ظ :(ظ) 


سیال و همکاران (2015 ,.21 66 81721) در پاکستان با استفاده از 
شاخص ]1۷ حاصل از تصاویر دوره هشت ساله سنحنده + ۳1۷۲ 
مقدار عملکرد برنج را در سال ۲۰۱۳ تخمین زدند. مدل رگرسیونی 
پیشنهادی آنان بین مقدار محصول برداشت شده و مقدار شاخص‌های 
1 در اوج دوره سبزینگی گیاه با ضریب تبیین معنی‌دار برای تخمین 
عملکرد به‌کار گرفته شد. کاربرد شاخص ]1۷ برای پیش‌بینی عملکرد 
برنج (2019 .2 )6 عمه ۷۷ :2019 ,.21 6۲ ۷20۳01011 )» رابطه معنی - 
دار شاخص ]1۷ با محتوای آب پوشش گیاهی گندم ,21 6 «دنک) 
(2013. همبستگی بالای بین تجمع نیتروژن برگ گندم و برنج با 
شاخص ۷۲ (2008 ,.21 66 2۷)» برآورد عملکرد وش پنبه با استفاده 
از شاخص 5۷۲ (2012 ,.21 61 00167762))» همبستگی بالای شاخص 
۲ با عملکرد روغن پالم (2013 ,.21 6 0۲270صاو12ظ) و کاربرد 
شاخص ]1۷ در بررسی رشد و عملکرد نیشکر 50661»۳71) 
(م 076171077 با ضریب تبیین معنی‌دار (2005 ,.لة 61 150065) 
گزارش شده است. رابطه رگرسیونی عملکرد ثبت شده و پیش‌بینی شده 
با کمک شاخص ]1۷ در ۸۰ نقطه آزمون در سال ۱۳۹۵ و ۸۶ نقطه 
آزمون در سال ۱۳۹۶ نیز رابطه معنی‌دار بین شاخص 1۷1 و عملکرد 
دانه برنج را نشان داد (شکل ۲). شکل ۳ عملکرد واقعی برآورد شده 
برای استان گیلان را بر اساس مدل منتخب نشان می‌دهد. در سال 
۵ بر اساس شاخص ۲۷۲ عملکرد واقعی در اراضی زیر کشت برنج 
دارای دامنه ۲۷۵۵ تا ۴۶۵۰ کیلوگرم بر هکتار بوده است که بیشترین 


آن در شهرستان آستارا (۴۱۲۵ کیلوگرم در هکتار) و کمترین آن در 


شهرستان شفت (۲۹۹۱ کیلوگرم در هکتار) مشاهده شد (شکل ۲. 
همچنین در سال ۱۳۹۶ بر اساس مدل منتخب. عملکرد واقعی دارای 
دامنه ۲۲۷۰ تا ۴۷۱۴ کیلوگرم در هکتار بوده است که بیش‌ترین مقدار 
به شهرستان تالش (۳۳۱۶ کیلوگرم در هکتارا و کمترین مقدار به 
شهرستان شفت (۴۲۶۰ کیلوگرم در هکتار) تعلق داشت (شکل ۲). 


بررسی عملکرد واقعی با استفاده از روش زمین‌آمار 
بررسی ساختار مکانی داده‌ها: به‌منظور بررسی ساختار مکانی 
داده‌های مورد استفاده ابتدا مناسب‌ترین مدل به داده‌ها برازش داده شد 
([شکل ۴) و سپس نسبت بین اثر قطعه‌الی به سقف مورد بررسی قرار 
گرفت. وقتی این نسیت کمتر از ۲۵ درصد باشد متفیر مورد نظر ساختار 
مکانی قوی دارد اگر بین ۲۵ تا ۷۵ درصد باشد داده‌ها ساختار مکانی 
متوسط دارند و زمانی که بیشتر از ۷۵ درصد باشد ساختار مکانی ضعیفی 
دارند (1994 ,بلح )6 2۳002706112)). نسبت اثر قطعه‌ای به سقف در 
متغیر عملکرد دانه ۶٩(‏ و ۲۷ به‌ترتیب در سال‌های ۱۳۹۵ و ۱۳۹۶) 
بیان گر همبستگی مکانی متوسط و مناسب این متفیر است (جدول ۶ 
مقدار کم‌تر نسبت اثر قطعه‌ای به سقف و دامنه تأثیر بالاتر در سال 
۶ نسبت به سال ۱۳۹۵ نشان‌دهنده همبستگی مکانی بهتر داده‌ها 


در سال ۱۳۹۶ نشان‌دهنده پایین تر بودن خطاهای اندازه‌گیری در منابع 
تغییر مکانی در فواصل کوچک‌تر نسبت به نقاط سال ۱۳۹۵ است 
(جدول ۶؛ شکل ۴). 


۸ نشریه بوم شناسی کشاورزی, جلد ۱۴ شماره۳ . پاییز ۱۴۰۱ 


( "فطع 
26( 
262 
239 
270135 
7-2-6 
26299875 
27003 
2۰-01 
7-2-4 
2-0 
7-2-1 
2-194 
223-164 
27 
2-3 
7-2-76 
2206810 
7-2-2-0 
262-103 
221-2683 
2-2-8 
7-2-11 
7-22-38 
22-327 


جدول ۶- روابط خطی بین شاخص‌های گیاهی و عملکرد در سال ۱۳۹۵ 
6 ۱ ۲۱۱۵۱0 مضه مععن4ص ممتاهاع ۲6۵ 06 جهمع660 دمرتطاددم) ماع متا م1 -4 ماه 1 


جدر میانکین ری 
فرس ‏ 9 0 
۰ «مناحمنصهاه 
572 006 
7171 005 
572 0073 
572 006 
570 000125 
32 006 
572 006 
5.9 0023 
598 036 
598 035 
596 +034 
598 036 
3172 004 
5.9 037 
598 036 
5.390 97 0 
5.78 0019 
582 020 
575 **0۳12 
587 018 
58 003 
538 018 
6 0۳17 
536 018 


معادله 
۱۱۹۹ 


1098.5 +- 1 
1135.17 + 6 
14.953 + 6 
2001.1: + 7 
1098.5: + 1 
962.43 + 2 
1652.0: + 8 
13.044 + 5 
2655.77 +- 3 
3142.4 + 5 
25.864 + 2 
4765.27 - 4 
-17.775 + 3 
-2244.4( + 4 
3995.77 + 8 
34.342 + 6 
1720.5 + 5 
2516.77 + 4 
15.369: + 6 
2774.4 + 6 
-14.447 + 8 
29.701 + 8 
2323.4 + 6 
-1306.9( + 3 


1۷ 
۷ 
(1-9-۷ 


تاریخ تصویر 
مص)260 1۳۵220 
1916 


۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
۹ خرداد (18 10106) 
٩‏ مرداد (9 2۷6۲ع۸) 
مرداد (9 ۷9۲ع۸) 
٩‏ مرداد (9 ۷۲ع۸) 
٩‏ مرداد (9 ۸2۷۲) 
٩‏ مرداد (9 ۷9۲ع۸) 
٩‏ مرداد (9 ۸2۷6۲) 
٩‏ مرداد (9 ۷9۲ع۸) 
٩‏ مرداد (9 ۸2۷۲) 
۱ مرداد (21 5۲تاعتاض) 
۱ مرداد (21 5۲تاعباض) 
۱ مرداد (21 5۲تاعتاض) 
۱ مرداد (21 5۲اعباض) 
۱ مرداد (21 5۲اعتاض) 
۱ مرداد (21 5۲تاعباض) 
۱ مرداد (21 5۲اعباض) 
۱ مرداد (21 5۲تاعتاض) 


#و :25 به ترتیب معنی دار در سطح احتمال ۱ و ۵ و غیر معنی دار. 
۰ بر 21) لصم زو- مه 4ص واتانطاه‌طامتس 56۵ 4ص 19۵ 0۶ 16۷۵1 مطا 26 اصم‌تگتصعزه رفظ 200 بر 


تعداد 
نقاط 

۷0۰ ۶ 

ت۱۱ 
309 
309 
309 
309 
3209 
309 
309 
309 
309 
309 
309 
309 
309 
309 
بزازه 
309 
بزازه 
309 
309 
بزازه 
309 
309 
309 
بزازه 


ماهواره 
۱۳ 


8 200926[ 
لندست۸ 


8 200626[ 
لندست۸ 


8 2800521[ 
لندست۸ 


سال 
22۳ 


2016 


جذر میانگین 
مربعات خطا 
5۳( 
( 02.02 
1۱99-6 
199.66 
1۱93938453 
192976 
20-8 
1990668 
1۱32-51 
20032 
200.46 
200-31 
20038 
7-2-9 
1914 
200-1 
20035 
19626 
1۱-91 
19293-30۳69 
1۱-1027 
1939223 
19314 
19139944 
1۱2-61 
19-1347 
1۱-4121 
1۱93-952 
1۱93-923 
1۱-11 
1۱93391917 
1۱-2-67 
1۱99-46 
199.066 


جذر میانگین 
۵( 
(صعه‌ه۵۵ 
3295 


4.00 
3293 
3296 
32.92 
3296 
3296 
32.92 
32.92 
32.92 
32.92 
32.92 
2917 
32.92 
32.92 
32399 
105 
105 
112 
6 
32.99 
107 
105 
1-13 
1-18 
119 
1-18 
119 
2394 
119 
1-18 
1-13 


۲ ات00۵ 
«متاحصنصه/ع0 


)00۳04** 
)006** 
0)0023** 
)00۳0۹** 
00015 
)005** 
)005** 
00015 
00015 
00015 
00015 
00015 
)003** 
00015 
00015 
)003** 
00215۴ 
)0015** 
)0(218** 
00215۴ 
0010 
)0(16*#* 
0015۴ 
)0(16*#* 
)032** 
031+ 
)032** 
)032** 
)009** 
)33** 
)032** 
0)021** 


۰ ,20 رطع 1و -صمه مه اتازمادطامنم 59۵ صع 160 0۶ 1۵7۵1 عطا 24 اصهم‌تگنصعنو :ده 4ص ۲ 


معادله 
۱ 


313.737 + 9 
462.64 + 5 
9.6155 + 2 
637.877 + 89 
-3.9427 + 6 
-361.27 +4 90 
474.06: + 09 
-7 0839 + 2 
36.548 + 9 
-5.6567 +4 2 
1.01017 + 7 
65.075: + 9 
7.57647 + 2 
-30.2665 + 9 
148.717 + 5 
3.52765 + 7 
584.22 + 5 
777.99 + 4 
2.3753 + 1 
1155.47 + 9 
3.23487 + 4 
-622.187 + 9 
878.317 + 3 
7.2577 + 7 
1549.3 + 1 
1951.47 + 0 
16.667 + 2 
2731.9 + 8 
-15.2757 + 5 
- 1254.37 + 5 
2332.4 + 2 
22.443 + 1 
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06 


۱۱۸ 
۱۱۱ 
1۷ 
1۷ 
۳۷ 
1۷ 
۷ 
۳-3-2 
۱1(۷ 
۱۱۱ 
1۷ 
1۷ 
۳۷ 
1۷ 
۷ 
۳-3-2 
۱1(۷ 
۱۱۱ 
1۷ 
1۷ 
۳۷ 
1۷ 
۷ 
۳3-2 
۱۱(۷ 
۱۱۱ 
1۷ 
1۷ 
۳۷ 
1۷ 
۷ 
۳-3-2 


جدول ۵ روابط خطی بین شاخص‌های گیاهی 9 عملکرد در سال ۱۳۹۹ 
7 1 ۷۱۱۵۱0 0صه مععنهص هناهام عع۲ 6 م60 دمرتطافدم) ماع بمم‌ناً معط 1 -5 ما20 1 


تاریخ تصویر 
6 م260۳ 1012826 


)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
)1006 30( تیر‎ ٩ 
شهریور (13 06۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
شهریور (13 06۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
شهریور (13 006۲ع)960)‎ ۲ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
)11017 23( ۱مرداد‎ 
مرداد (8 )9تاعبنض)‎ ۷ 
مرداد (8 9۲تاعباض)‎ ۷ 
مرداد (8 )9تاعبتض)‎ ۷ 
مرداد (8 )9تاعباض)‎ ۷ 
مرداد (8 )9تاعبتض)‎ ۷ 
مرداد (8 9۲تاعبنض)‎ ۷ 
مرداد (8 )9تاعباض)‎ ۷ 
مرداد (8 )9تاعبنض)‎ ۷ 


#و :06 به ترتیب معنی دار در سطح احتمال ۱ و ۵ و غیر معنی دار. 


تعداد نقاط 

۱ 

عاصنم ]0 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 
350 


ماهواره 
۱ 


ماعصتاجهد 


سنتینل ۲۳ 


ماعصتاصهو 


سنتینل ۲۳ 


8 200521[ 
لندست ۸ 


8 200521[ 
لندست ۸ 


2۸۹ 


سال 
91۳ 


2017 


2017 


2017 


2017 
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17( 26 
5000 
10718-471 ۳۹ مس 
* 2 7۲ ۳ ۶ 
جع موم ظ گِ 
ظه ظد 
‌ ت 
۵ 9 
:۳ :۲ 
4000 کل ۲ 
5 9 
0 ب> 3 
2000 
5000 00 4000 350 200 40 200 3700 3200 
(۲۵.21) ۷210 24) :۳:۵ (8.521) ۷۱10 ۳:2۵:20 


شکل ۲- رابطه رگرسیونی بین عملکرد دانه ثبت شده و پیش‌بینی ده الف) با کمک شاخص ]13۷ در ۸۰ نقطه آزمون در سال ۱۳۹۵ و ب) با 
کمک شاخص 13۷۲ در ۸۲ نقطه آزمون در سال ۱۳۹۲ 
0010/5 )69 80 2 ۱۵06۵ ]۲۴۲۷ ۵۶ ماع 6 6 ره ۲۱۵۱0 صقرم ۵۳۵016660 200 ۲۵۵۵۲0۵60 م۵۲۵۵ مرتصاعصمتها۲۵ «ملدوم عم - 2 ۲۱۵۰ 
7 ۱ 001069 691] 86 2 1۱۳06۲ ۲۰۲۷۲ ۵۶ ماع عظ) ۳۷16 (ظ 20 2016 


000 00 3-0000( 3-0000 0000 000+ 3000 2-0000( 
ِ ِ ِ_ هگ ِ سح ِ_ 
2017 2016 
هِ 2 ِ ۰۵ (ع 
۳ دو ۳ ده 
3 3 چِ چ 
[ه ‌ِ ‌ ‌ 
تب تب تب ف‌ 
هِ هِ 
۳ 3 8 1 
3 و جع ۰ 
ك ه 3 3 
چ چ چ چ 
۳ 2 ۳ اافت 
ِ ِ ِِ ِِ 
اک ادا ادا لح 
چ ت چ 
۵ ماش 
ماو 
هِ هِ 
ِِ 5 5 5 
ّ 1 ِ 0 مس 21 
ِ 5 ِ 
چ چ چ چ 
255 
لت تس تلا ۰۰ص 
0 80 40 0 10 0۰ 80 0 40 0 10 0 
45000 000 350000 30000 450000 400000 350000 200000 


شکل ۳- نقشه عملکرد واقعی استان گیلان بر اساس مدل منتخب مبتنی بر شاخص گیاهی 1*۷۲ در سال‌های ۱۳۹۵ و ۱۳۹۲ 
7 2:0 2016 10 ۲ ۲۸۱۷ 66 چم 60و۵2 0۵مص ۵660و 6 چم 0مووه معصز ۵ص صهانیتی ۵۶ مقر ۵10 [۸۲۵ -3 ۲۱۵۰ 
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2( 6 0( 17 


0 ۱2) عممداوزل تعوجع ‏ م8 ۰ ۱ سس ۱۰107 (عاع/) عمجهاوای فعوه ‏ 8 ۰ م۵ماا سس 


۷۳105 ۷۳*103 
ِ ۵2 ۵ ۰ 32839 
0819 ۳ 289 
054 ی ی تس ۳ 
93 #6" ال ّ 096 


‌ 
7 5 29 1.152 0.576 0 9 4158 3327 2.495 1.653 0.832 0 
شکل >- نیم تغییرنماهای تجربی و مدل برازش داده به آن‌ها با روش-مدل کریجینگ معمولی 902016 برای عملکرد دانه در الف) سال ۱۳۹۵ و ب) 
سال ۱۳۹۲ 
صتجه 10۲ 0-006م)مصه ماطاهاز عصنعاه۱! وتحوتابن عط) زوا معط 6۵ 110660 061مصه م6 20 قصصه و مزنه نهد ه)صمص۵۵ 1۲ -4 ,و۲ 
7 (0 200 2016 (2 صا ۷1۵10 

# خط آبی مدل نیم‌تغییرنما برازش داده شده نقاطء نقاط قرمز نقاط جمح‌آوری شده علامت‌های مثبت معادل میانگین نقاط نزدیک به هم‌جهت جلوگیری از تداخل نتایج است. 
مصعتو ۷۵تاتووم ما فصه رعاصتمح 6۵116060 مطا کجمومتم1۵ عامل ۲۵0 مط رقاصتمم 060 عمط ۵۶ معط ماماهتیه۲-تصعو مطا ۵۶ مصناً میتاها 1۵ بو 


من فالناوع قط) غصمبعنج ما قطان طمعع ما مومام فاصتمم مفط ۵۶ مهههع27 عطا ما لقناوم 


جدول "- مولفه‌های نیم تغییرنماهای مورد استفاده در مطالعه برای بررسی ساختار مکانی داده‌های عملکرد برنج در سال‌های ۱۳۹۵ و ۱۳۹۲ 
0 2016 1۱ 0212 ۲۱۵10 م۱6 ۵ هساو آهتادصه م) )و1۳۲ ۵) وی م6 جز 60و مهو هدعو ۵ دم2۵ ۴2۲ -6 12016 


2017 
اثر قطعه‌ای به سقف دامنه تأثیر سقف اثر قطعه‌ای مدل سال 
ایهم 5 ما ۲۳۲6۵ امعع۱۱ 10926 الا 61 ۱۱۵260( [۷۲006 1۳( 
69 3350 329261 5644 ثابت )٩)2016(‏ 2016 
27 20249 3239803 ثابت )٩)2016(‏ 2017 


روابط خطی به‌دست آمده» ضریب تبیین و ضریب تغییرات مرتبط 
با ۳۲ روش-مدل درون‌یابی مورد مطالعه و عملکرد در سال ۱۳۹۵ در 
جدول ۷ نشان داده شد. در سال ۱۳۹۵ بر اساس کم‌ترین جذر ميانگین 
مربعات خطا (۳۰۱/۵۵) و ضریب تبیین معنی‌دار روش-مدل کریجینگ 
ساده -1-565561 بهترین روش-مدل در درون‌یابی عملکرد دانه بوده 
امّا به‌دلیل عدم مشاهده روند مشخص از روش-مدل کریجینگ 
معمولی-62016 جهت درون‌یابی عملکرد استفاده شد (جدول ۷. 
همچنین در سال ۱۳۹۶ بر اساس کم‌ترین جذر میانگین مربعات خطا 
(۳۰۱/۵۵) و ضریب تبیین معنی‌داره روش-مدل کریجینگ معمولی- 
6 بیش‌ترین دقت را در درون‌یابی عملکرد دانه برنج داشته است 
(جدول ۸). 

نتایج این پژوهش برتری کریجینگ معمولی نسبت به سایر روش - 
های دزون‌بابی شامل کریچینکت جهانی و کریچینک ساده را نشان داد 
در تحقیقات دیگر نیز کاربرد کریجینک معمولی در درون‌یابی عملکرد 
برنج و خصوصیات خاک در جنوب چین (2014 ,.21 61 ا1) و کاربرد 


مطلوب کریجینگ معمولی در برآورد نوسانات عملکرد مزارع گندم 
(2016 .له اه تامآمطه]۷ تتنعد) گزارش شده است. این در حالی 
است که اژیرابی و همکاران (2019 ,۵1 اه نطهتنطتش) در برآورد 
عملکرد گندم در مزرعه نمونه ارتش آققلاء برتری مدل‌های کریجینگ 
جهانی نسبت به کریجینگ معمولی را گزارش نمودند. تسفاهونگن و 
همکاران (2011 ,له 6۲ 1651]200620) علت برتری کریجینگ 
معمولی نسبت به سایر روش‌های درون‌یابی را ارائه بهترین نتایج جهت 
درون‌یابی فواصل نمونه‌گیری نشده و به حداقل رساندن اثر نقاط پرت 
روی نتایج درون‌یابی گزارش نمودند. رابطه رگرسیونی عملکرد ثبت 
شده و پیش‌بینی شده با مدل ثابت در روش کریجینگ معمولی ۸۰ 
نقطه تست در سال ۱۳۹۵ و ۸۶ نقطه تست در سال ۱۳۹۶ نیز رابطه 
معنی‌دار بین مقادیر درون‌یابی شده در روش کریجینگ و عملکرد دانه 
برنج را نشان داد (شکل ۵). شکل ۶ عملکرد واقعی برآورد شده برای 
استان گیلان را بر اساس روش-مدل منتخب نشان می‌دهد. 


۲ نشریه بوم شناسی کشاورزی. جلد ۰۱۴ شماره ۳ . پاییز ۱۳۰۱ 


جدول ۷- نتایج ارزیابی ۳۳ روشمدل مورد استفاده برای درون‌یابی متغیر عملکرد دانه برنج در استان گیلان [ کریجینگ معمولی (016)» ساده 
(51) و جهانی 1/160)] بر اساس جذر میانگین مربعات خطا (180151) در سال ۱۳۹۵ 
۸۵ 6۱۷۱۱2 رز ماصاهزیه۱ ۷۱۵۱0 صنواع م۱ ۵۶ صمتهاهصه)صا 10 0عوتا و0 00۱-0 مه ۵۶33 ممت)مادی ۵ داوم -7 12016 
6 1۳ ۷۲۹۲ ماه 50۵۵۵ همه )۵۵ 606 صم موه (عاتنا) اهتنا 2۸0 باه رصن را متفصتل 0۵] 


جدد یکت .نیراد زقایل 
مرعات ها ۶ باحصا 
۱ 1 ۳0( 
32196 320 
323.3 320 
29542 320 
320991 320 
322447 320 
32113 320 
299.93 320 
29978 320 
20034 320 
320112 320 
320117 320 
30242 320 
320105 320 
3202.5 320 
29983 320 
2384 320 
2-38 320 
320432 320 
20298 320 
32028 320 
2014 320 
205 320 
294106 320 
25 320 
29017 320 
20961 320 
297.59 320 


جذر میانکین 


۳ 
(۵0ع-۵۵) 


214 
32.17 


290 
325 
3213 


214 


2.70 
209 


209 


208 
201 
238 
25 
237 


2.12 


242 


244 


۶۴ اصمن‌تکله) 
ما حصنصصه6ع0 


ش۱0۳0 


0 0351+ 


0042** 


0۳0 


0055** 


0057+ 


)0.031** 


0030+ 


0030+ 


0030+ 


0030+ 


0095 


012** 


«جهوم.0 


011 


ب+«+03096۹ 


0095 


0082+ 


0091+ 


0092+ 


0093** 


0094** 


0045 


سز9 0۳0 


)0044** 


005 


0030+ 


معادله 
۱۹ 


+ * 00980 
201.723 
+ * 00943 
3608441 
+ * 00501 
2364241 
+ر* مهو 
720-7 


0.1083 * + 
(3-6 
0.0980 * + 
2-3 
0.0069 * + 
32-59 


0.0058 * + 
2-3-14 


0.0058 * + 
2-2-2-4 
0.0080 * + 
2498 
0.0081 * + 
2-30 
0.0311 * + 
29068 
0.0903 * + 
265.910 
0.0280 * 4+ 
32152 
0.0525 * + 
38029399 


0.0290 * ٩4+ 
2-2-36 
0.0311 * + 
(0 
0.0288 * + 
(3-1 


0.0318 * + 
(2-2 3 


+ * 0.0323 
2-62 
+ * 0.0331 
2-118 
+ * 0.0345 
2-2-39 
+ * 0.0432 
2-8( 
+ * 0.0942 
2092470 
+ ۶ * 0.0438 
2-34" 
+ بر * 0.0893 
2-2237 


+ بر * 0.0084 
3-01( 


مدل 
۷00۵ 
ثابت (812016) 
([556عظ1-۳) 


)1۳6556[( 


اثر روزنه‌ای ۳۲016) 
(۳1]60۱ 


درجه دو منطقی ۲0110021) 


(410ه۲او 
گوسی (02055120)) 
نمایی (2۵0۴6121ظ) 


کروی پنج وجهی 
([162مصامکه]جعظ۳) 


کروی چهاروجهی 
(لهءتعطاموهتاع) 


کروی )٩0۳6:6۵1(‏ 
دایره‌ای (ت12ا0170)) 
ثابت (512016) 
([56وع1-۴) 


(اعوعع1-۳۴) 


اثر روزنه‌ای ۳۲016) 
(۳1]60۱ 


درجه دو منطقی ۲010021) 


(410ه۲او 
گوسی (02055120)) 
نمایی (لهتاجهعمحظ) 
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کروی چهاروجهی 
(لهءتعطاموهتاع) 


کروی )٩0۳6:6۵1(‏ 
دایره‌ای (تق1ا0170)) 
ثابت (812016) 
([6556ظ-1) 


)۲۳6556[( 


اثر روزنه‌ای ۳۲016) 
(۳]]60۱ 


درجه دو منطقی ۲010021) 


(410ه۲او 


روش درون‌یابی 
ممتاها0نه)1 
00( 


01 
۵1 
۵1 


۵1 


سال 


21۳ 


2016 


2016 


219۹44 


29038 


29789 


2.054 


22934 


22904 


320 


320 


320 


320 


320 


320 


214 


215 


25 


25 


25 


25 
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0.0979 * + 


اد چِ 9 
057 0۳ 1247 2609 گوسی (02055120)) 
ید + * 0.0525 ۲ ۲ 
0033 2795640 نمایی (لهتاج6عظ) 
003۹0۷ + * 0.0079 کروی پنج وجهی 
32-1 (1621عصامکه]جعظ۳) 
0030۷ +۶ * 0.0079 کروی چهاروجهی 
زازز رز (لهءتعصاموه0ع1) 
"۳ + * 0.0079 ۲ 
0029 5018 2957 کروی )٩0۳61621(‏ 
"۳ + * 0.0080 ۰ 
0۳29 4544 3957 دایره‌ای (ع۳0012)) 


#و :06 به ترتیب معنی دار در سطح احتمال ۱ و ۵ و غیر معنی دار. 


۰( راهم زانطع زوجم مج زاتانطاهطامج 5۵ 4صح 16۵ 0۶ 16۷۵1 عط) )2 )صهم‌نصعله :6ط 0ه2 


۳ 


2016 


جدول ۸- نتایج ارزیابی ۳۳ روش-مدل مورد استفاده برای درون‌یابی متغیر عملکرد دانه برنج در استان گیلان [ کریجینگ معمولی (016)» ساده 


(516) و جهانی (116)] بر اساس جذر میانگین مربعات خطا (:1۷151) در سال ۱۳۹۲ 


۸۵ 6۱۱۱2۵۳ صز ع۱صاهیه۲ ۲۱۵۱0 تج رمع ۵ حمتاهاهمعماصز دم 0موت ولم)ع-00 مه ۵۶33 ممتامادی 0 ددع1 -8 م1201 


7 1۳ ۱۷۲۹۲ ۵۵۲ 50۵۵۵ مهم )۳۵۵ 606 دی 0عفهه زا اهتنا مضه رام مامرصنه را فطل 0۵۳ 
جذر میانگین 
مربعات خطا 
۱/۵ 
(اجه۵ع) 


جذر میانگین 

مربعات خطا 
۳ 
(تمط.ع) 


2271 
22.6 
22-1 
2294 
226217 


28 


21982 


2-4 


2-1-6 


21-68 
21917 
228 
2305 


2281 


2137 


22023 


25 


24 


229.06 


تعداد نقاط 


۱۱0۵۲ 
رناا ۱۱ 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


340 


238 


2.40 


247 


241 


241 


241 


243 
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مقایسه رهیافت‌های مبتنی بر سنجش از دور و زمین‌آمار و 
برآورد عملکرد واقعی 

مقایسه روش‌های مبتنی بر زمین‌آمار و شاخص‌های گیاهی به‌دست 
آمده از تصاویر ماهواره‌ای در برآورد عملکرد واقعی دانه برنج نشان داد 
که بر اساس جذر میانگین مربعات خطا شبیه‌سازی و ضریب تبیین 
به‌دست آمده بین مقادیر واقعی و شبیه‌سازی شده عملکرده شاخص 
1 حاصل از تصاویر ماهواره‌ای ۱٩‏ مرداد در سال ۱۳۹۵ و ۱۷ مرداد 
در سال ۱۳۹۶ به بهترین شکل مقادیر عملکرد برنج محدوده مورد 
مطالعه را برآورد کرده است (جدول‌های ۴ ۵ ۷ و ۸). محمدی 
احمدمحمودی و همکاران ۲ نتامحصط۱۷2 ممصطم تمصصصعفطم]۷() 
(2015 ,.له در مقایسه روش زمین‌آمار و استفاده از داده‌های سنجش 
از دور به‌منظور پیش‌بینی عملکرد گندم در برخی از مراحل رشد در 
استان گلستان» در ارزیابی نقشه‌های عملکرد پیش‌بینی شده, برتری 
قاطع روش‌های مبتنی بر سنجش از دور را نسبت به زمین‌آمار نشان 
دادند. تایح این پژوهش بین‌گر قابلیت تصاویرماهوار‌ای در پیش‌بینی 
عملکرد گندم در اواخر پنجه‌زنی با ۷۱۵ کیلوگرم در هکتار خطای برآورد 


بود. تفاوت به‌دست آمده در دقت دو روش زمین‌آمار و مدل-روش‌های 
مبتنی بر سنجش از دور می‌تواند به‌دلیل تفاوت اساس دو روش در 
برآورد متغیرها باشد. کریجینگ یک روش زمین‌آمار برای درون‌یابی 
داده‌ها بر اساس واریانس مکانی است که به هر دو معیار فاصله و میزان 
تغییرات در زمان پیش‌بینی مقدار مجهول توجه می‌کند ۲06۲2۷1) 
(2014 .21 ا6. این در حالی است که در روش‌های مبتنی بر سنجش 
از دور عموماً شاخص‌های گیاهی برای ارزیابی یک یا چند متفیر خاصض 
بسیاری از شاخص‌های گیاهی بر مبنای بازتاب در محدوده نور قرمز و 
مادون‌قرمز نزدیک قرار داشته و ترکیبی جبری از این دو باند هستند 
(1995 ,صعطل > »۲۱۷102 :2008 رل 6 ت22ظ). بر اساس 
روش مبتنی بر سنجش از دور داده‌های باندهای طیفی مختلف می- 
توانند درصد پوشش متغیرهایی مانند گیاهان» زیست‌توده و شاخص 
سطح برگ را در هر واحد سلول تصاویر و با دقت مطلوب پیش‌بینی 
نمایند (2008 ملع 6 820216060[20). 


در اين پژوهش, از داده‌های نمونه گیری شده توسط ناظرین برنج 
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در سال‌های ۱۳۹۵ و ۱۳۹۶ استفاده شد که پراکنش نقاط نمونه- 
گیری‌شده در شکل ۱ نشان داده شده است. یکی از دلایل کاهش دقت 
روش کریجینگ در پیش‌بینی عملکرد دانه در محدوده مورد مطالعه 
می‌تواند تعداد نقاط نمونه ناکافی و بر اساس پراکنش نامناسب باشد که 
سبب می‌شود تا روش کریجینگ تغییرات عملکرد اراضی را با دقت 
کم‌تری پیش‌بینی نماید. ون‌گرونیگن و همکاران 00601867 (۷۵) 
(1999 ,.21 66 گزارش کردند که الگوی مکانی نمونه‌گیری در روش 
زمین‌آمار می‌تواند یکی از دلایل مهم در به حداقلرسانی واریانس روش 
درون‌یابی کریجینگ باشد؛ بنابراین می‌توان نتیجه گرفت که اجرای 
تموفه گیری دز سح اسان با اسفاده از الگرتن مضخص با بدکارگیرش 
روش‌هایی مانند اجرای الگوی عصحعت1 در ۸0۲5 در محدوده 
کاربری تخصیص یافته به کشت برنج» می‌تواند سبب افزایش دقت 


درون‌یابی با روش زمین‌آمار گردد. 


نتيجه گیری 

رهیافت‌های زمین‌آمار و مبتنی بر سنجش از دور از روش‌های 
پرکاربرد در پیش‌بینی عملکرد در یک محدوده می‌باشند که می‌توانند 
تارج کاربردی و مفیدی از موجودیت عملکرد در یک منطقه جهت انا 
سیاست گذاری کلان در تولید را فراهم نمایند. نتایج نشان داد که پایش 
عملکرد در سطح استان گیلان» می‌تواند اطلاعات مفیدی در جهت 
پیش‌بینی میزان عملکرد در کل استان را فراهم نماید. در مقایسه روش 
زمین‌آمار و روش-مدل‌های مبتنی بر سنجش از دور در پیش‌بینی 
عملکرد مزارع مشاهده شد که روش مبتنی بر سنجش از دور بر اساس 
نمونه گیری‌های انحام شده دقت بالاتری داشته است. روش زمین‌آمار 
و سنجش از دور دو روش متفاوت در برآورد یک متفیر در سطح بر 
تین اسان آتداز گیری اه قطانم تاک ووژن امه مار وی 
تغییرات متغیر را بر اساس معادلات ریاضی پیش‌بینی می‌کند. این روند 
می‌تواند تحت تأثیر عوامل مختلفی از جمله تغییر آب و هوا؛ خصوصیات 
اک ما اند پل و اه اس ان کي اش ار گر 
همچنین استفاده از روش زمین‌آمار جهت اجرای یک درون‌یابی مطلوب 


نیازمند نمونه گیری بسیار منظم می‌باشد که اغلب با توجه به کمبود 
امکانات امکان‌پذیر نمی‌باشد. روش سنجش از دور را می‌توان از اين 
حیث مورد توجه قرار داد. زیرا روش سنجش از دور با توجه بر مبتنی بر 
انعکاس امواج الکترومغناطیس بودن» دارای مزایایی از جمله امکان 
برآورد متفیرها در مناطق با دسترسی دشوار و عدم تأثیرپذیری از 
تغیرات محیطی می‌باشد. مزیت این دو رویه در مطالعات آگرواکولوژیک 
مقیاس وسیع فراهم می‌کنند و در نتبحه» استفاده از این اطلاعات را در 
مدل‌های مکان‌مند تصمیم‌سازی" ممکن می‌سازند. این روش‌ها به‌ویژه 
در شرایط که هدف مقایسه نسبت‌ها و يا استفاده از شاخص‌های نسبی 
باشد در مقایسه با شرایطی که هدف تعیین مطلق مقادیر یک کمیت 
باشد بسیار کارآمد هستند و می‌تواننند در حوزه مطالعات آگرواکولوژیک 
راهگشا باشند. 

پیش‌بینی پیش از موعد برداشت بر اساس اطلاعات مستخرج از 
این پژوهش نشان می‌دهد که استفاده از تصویر /تصاویر به‌دست آمده 
در ماه‌های خرداد و مرداد و سنجنده آآ۵ ماهواره لندست-۸ در استان 
گیلان می‌تواند به‌عنوان مبنای پیش‌بینی‌های عملکرد اين گیاه در 
سال‌های آتی مدنظر قرار گیرد. چنین مطالعاتی در استان گیلان, با 
توجه به سهم استان گیلان در تأمین برنج کشور نقش موثری در 
تصمیم‌گیری‌های مدیریتی در خصوص عرضه و تقاضای برنج در سطوح 
کلان خواهد داشت. همچنین نتایج نشان داد که نمونه‌گیری انحام شده 
در محدوده مورد مطالعه بر اساس الگوی مکانی منظم (سیستماتیک) 
موجب افزایش دقت روش‌های زمین‌آمار در تخمین عملکرد منطقه‌ای 
خواهد شد. نتایج این طرح می‌تواند اطلاعات پایه مناسبی برای مطالعات 
دیگر از جمله بررسی میزان خلاء عملکرد. دلایل خلاء عملکرد بررسی 
ارتباط تناسب اراضی با میزان عملکرد به‌دست آمده و پیش‌بینی و برآورد 
عملکرد در بازه زمانی مختلف را فراهم نماید. 
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15001-11۰ :58 داز تمد گه مهم معصمهگ لر60 .ولر۵و 10۳/۵ جع صا فمتاتممم۲0م 011و ۵۶ ۷2۲1۵01110 ۲۱۵۱0-90216 


)05://001.012/ 10.2 136/598521 76 

۵ 200 ۸۵۸۵۲۴ 08 0611۷60 »1۳06 2762 کعع1 ۵۶ ممت)هاتصطزووظ .2008 و۷۲ و۲۵21 0ص2 و.۲ ما۷2 ,.ت) ,مطتله)هو رما ,12606 
1395-7۰ :(۱1204 اصمصممتصظ گ۵ . مصتفجعو. مامصصمگ . ,م۷ ۵2 قح ۵ 20۵061 ۱۷۷۳۲۵۵۲-فکنل) مفصاً عاق 
3.-005://001.0182/10.1016/.1562-0 

۰ ۷6۵602010 ۱621-10112760 صرح ۲60 ۱21۲0۷۷-0200 لهج 0۲۵20-۵20 ۵۶ ط0وتتهم‌صصمن ,1995 ,.2 رصعط 220 و.نز.) ,۲۱۷1026 
005://001.0۲8/10.1016/0034-4257)95(00132-6 .38-48 :(54)1 اصمحصمم ظظ ۵۲ عصتعمهو امعم 

,اج 0۲۵/1۵۵ ۱۱۵://۱۷۱۸۷۱۸۵۵۵ظ حص 0عانو) فصمتنجل 0عانصرنا عطا ۶ه صمتاهمنصهع:ن0 احتتطانه‌تنهه م2 ۳۲۵0۵۵ .2018 ۲۸۵۰ 
1۳1۳2020 

(صهلونعظ ظ1) .۵ 257 .عمناهعناهانام همه .عصتعوعد م)محصفک ۵۶ م1فظ .2006 ,.۷ ,۴622۵ 20 ,.ظ رتصع]۳2 

۷ 00122 .11۳6516 ۷۲.۹6۰ .نام صهععم0 0۶ و۱018 )2عظ۱« م1 بواتوه۷ 1۵01 ۷۷۵۵۵ عصتاهعناوه 1۳۷ 2011۰ .۱۷۲ رتصه‌فهطت) 
(«تحصصصصاگ طوتلمطظ طاذه صفتویعظ م1) .۵ 114 بصعت رومعتامومک لهعبهاظ مه معمصمتمه آمتمانه‌تتعض ۵۴ 

0صضص عصتمیه فص صمتاهءلوعداه صهام‌عمه: هه بعبامن صقن .2006 ریکایه ۱۵۹2۷1 مضه وی دنه و رتاو 
صهلوه۳ 1) .265-277 :(1303 ط0ت2عوم )۲عوع1۳ مه 0صهامعصه .0عطوم) 2۵ 2000 (رجمتتاد وفعت) (1۳) معفصطا ماتالملهه 
(«تححصحصصت 5 طوتاعصظ ۷16۸ 

)۱2 ۷/10 وعع41ص1 عمصفام1۵116 ۹0۵6۵021 0۴ ۸۵۹0012108 .2012 ,.) بعطه ۷۷ 200 ریکایک و101 ریک بقمات0اا ب۱۷ ر2عتمتاتان 
2 135/0100301201 10.2 /۳۲05://01,:0۴۲8 .849-857 :(52)2 501668 وم .قم۲ا۵ع صهامدا ا ۷1۵1۵ )1 هه 

,295-9 :25 اصمصصممصظ ۵۶ . عقوم مامصفگ ۰( 8۸۱۷) 46ص ممتاهاععع۷ )۵109و ۸ .1988 بیکایظه هام۲ 
-005://001.01۲2/10.1016/0034-41257)88(90106 

امءتافتامادهعع ۵۶ ممتامتاه۲۲۷ .2012 ر.و بلتطعع20 20 رب متفقاقطه ب.ظ رتفلطهک رم متطقاده92۳۷ لامعمصصطه 1 رب متصه‌مفکز 
٩۵۱۱ ۰‏ ۷۷۵۱۵۲ ۵۲۵۷۱۵۵۵۰ صهاوع6۵1 ۵۶ وصه1 لمبانهتتعه ما فاصم نهمممهه ۵۶ عصتمم2 مه عصتاحصصتاوه 10۲ و0‌مطامعظ 
(«تفحصصصری و فتاعصظ اه نهتوم۳ م1) .201-218 :(22)1 

۶۲ ۱6 مصتصصعهاع0 ما 08مطاممظر لوملافتاماوممع ۵ ممفتنوم‌جوم) ,2014 ون مصهتام۱۷]21 مه ,۷ متصصقتصهعا۲ ویک 105172۷ 
0 ۶ امصتام؟ صمتصهت! بصفع تمصع ما مملاباطتمعنل ومته‌عمد عصنام0مصظ عم 0مطامصه صمتاهاهمتعاصا ماه متاقصصئاه‌مزه 
و20.1001.1.24763128:1393.3.8.9 ( تحص 5 فاعم طاذ۳ صفتویعظ م) .55-68 :(3)8 «0102عظ 

0212۵۱6 00و ]۵۵2 0۴ ۳۴۵016۷۵1 .2013 ویک رععا 200 .تا معصلال وگ ,۲۱۵198 میک مطفتا‌صاظ ب.1 بط20160! ول ملک 
808-۰ . :(1104 . اما عصتعصمد... عامصصف ...20 عمصقاعوهعی) .. بط .. .عععت4ص1 . صمتاماموه۷7 ۰ ۲۵02۲ 
109/10195 10.1 /۲05://001.012 

بط .ماه عه4ه۲ ملتامصتتهاهن عصتوید مسافتمجظ مد متقصصتادع ما طمهمتمم2 وع1تعو-عصصتا ۸ .2009 و.ط با2 ص۷2 220 ,.۲ متطنک 
4 0۳۴8 109/1 ۳1)05://001.0۲8/10.1 .2519-2527 :47 عطاومهو م)محصفمک مه 0605616866 و0 فصمت1۲28826 1 

۵۵00۵ )۱/1 عظ2۵081)01 ممتک ,2013 و.ت) بطاهتظ 200 ویک یحع .۷ ممهل۱۷ و۷ ,۷۵20 و تال وی با ,۷۷۰ ,150006 
:1 ,صمتاه0صز0ع0) صرح ممتاهبعودان) اتف لمنامم۸ گه امتمژ لمعمتامصعماص ماه له همع[ مرتامصتتهاهمآهنل هه 
6 2 1005://001.012/10.1016/[.1282 .568-576 

۵ 1۱۱6۵ ۵۴ ات۷2۵ لهتلهوه وه اصمصروومعوه مامصصله ۸ .2014 ,.ت) رعوصولنا 2۳00 ,.م) راصا و.ظ مها و.ل تمه .۷۷ ت20 .2 بتا] 
.39-7 .۰ :235 قصصع۵0عی .. بقصتطی) . طانامو... صاً.. وف دص 0 .. فمتاتعمميم .. لقعتصعفط .لاصو .۰ 0)60عع۱عو . هه 
7 260061۳۵22 [/05://001.0۲82/10.1016) 


۸ نشریه بوم شناسی کشاورزی. جلد ۱۴ شماره۳. پاییز ۱۴۰۱ 


,56-4 :143 دعوم و۵ ۳۱6۱0 .ولمم مهع ۷۱۵10 ومع 10۲ 2 ماتلامنهه ۵۶ مویا م1 ,2013 ,.ظ.(1 مالع0ا0.] 
8 (/009://001.01۲82/10.1016 

۵۱0)0-۵۵60067 902016۵0 عصلمنا ممتاقحصتاوع ۲۱۵۵-۷1۵10 مل2مو-1ه9م1ع۴۵ .2019 ,.ل ,۴۱۵0 200 ریک رعمی ی بناج لا ,. 1.۷۷ ,1۷۲2 
51-7۰ :(40)1 عصلوجعگ ع)محصف ۵۶ هتم لمعمتامصعماص بقعرمک طانامق گم رنه موه ۸ تهج وآنا۱۷]0 مه متلقصتاه 
1 0 .005://001.012/10 

مج ممتام‌صنوجم ماه بلعد۷ وم ۵۴ طمتام۳۳۵0۲ ,2005 ویگا. ] رته16 ۱۷۱۵۷ 200 ولا بعحقله .لا ولاکم و2 انا و مان وک ,1۷۲0 
۲۵۵۱0 .صتق۱ظ مصنطن ات۱۱۵ عمط جم 0212 0ععصعو «امامصصعء عصتوته اع0مصه ط) مه ۸-0 2۷ اد «مصمکه‌تگله مو۲ا م۱۷۵۲ 
2 0.1016 1009://001.012/1 .301-322 :(183)2 عصتااع۱۷۲00 

عصتقصعو مامصصمع مه -لهمتافتاهاوممع ۵۴ ممولندم‌جط0) :2015 وم رتطا۸ 220 ر.ظ وتفلصفک ر.ظ رتتامصطه]۱۷ 0محصطمه تمصصهطم]۱۷ 
صمادعا00 ر0ع11 طممممصصهل :ولناه عفد ۸) فععماه عصاهع ۵۶ مصومو صا صمتام0۲601 ۷۱۵۱۵0 ۷/۵۵۵۲ ط1 فلمطامصظ 0202-02860 
5 ( ۱۱1۱۱۱۵2۲۷ مفتعصظ از صهتوع۲ 12) .51-76 :(802 ممتامنال۳۲۵ 00۵ 0۲0۵۷1868(۰ 

عصلوه ۷۱۵۱۵0 )7۱۵۵۵ ۵۴ ۷2/1201110۲ ۸۱6۱0 منز عصتافصتو۴ .2016 ۷ متصقطه! مضه بیکایه متکامط00مک وی ,۱2021124 تتلووهل 
حوناعص طاز۳ صمتویعظ مآ) , 329-345 :(803 وما0عع0ه۸ ۵۶ امراه[ .ععبمآنام‌توه مملفلمع؟م ما طم2(0102 مه :عمصاممه۷ لمتله0٩‏ 
2 05://01.0۲8/10.22067/122۷ظ (وتمصصت 

حچ عصلوتا وقع ۷1۵10 ۱۵۲276 101۲226 0۴ عمتلهیاله۲۷ .2019 ,یط رتنامم‌تقاطآ۱۵ 220 ریخ رتصمااهی و.ظ وتهلصهک وب مصهنمصتتا۳۵ 
صماوع00۱) طا مصلفدها عماج تناو :واه عفدن)) ‏ 0مطافجظ . اتم‌عمصضا ماتالمامعاع0مصه ممتاق‌اناصته ۵۵ 1۳]6812060 
 ( 9‏ 7 .253-268 :(65)2 م81 [ر50 ده چهمهمتو۸ ۵1 وعتطهت .(۲۵۷1868ظ 

ص موی ات۵ممنای مضه ردام رم عمتجم .2019 .لا رعش مضه 9.۸ مصللتاتفطقط: ریکلیک رقلاع0ع۲ ریگ متفصطف۴ 
۱۷۸0 :(5)1 تاه نت۸ 220 ۳00 عع0) .صتقطاه ایرفایی ااکزی ۱ ۱2۹/۸ 1۷۲2122 
2 ۱ 0 0۹/0( 

طاتحظ .نمی 20۳ همنسل مععصفط ما گنه ممتمع‌ناده ها 0صه مد ۵۶ ممنامی‌تاتوعهآل عامصنل مونعن-معع1۵0 .2017 ۲.۱ رلعزجم۴ 
6 ۳۲5۳۳۱ [/10.22059 /2 001۰0۲ //:۱۲)۵۵ .۰( تیاو افتاع ص طز صهتویع۴ ع1) ,419-439 :43 وعلووطظ م۵2 20 

(8۸۷) :106 عمناهاعع۷6 ٩0۵11-20[05]60‏ ص1 ما26 موز ۵11و ۱۵220۷6 عطلونا .20168 .۲ ,عطقم 220 ر.ت) بتا20 بلط مرطفگ 
.4139-5۰ :209 عصمصصمماظظ که عصلفصعگ ما)محصفک .وصماووجع تیه ما ممتتقصتاده ففعصصماها. فص 0صنامتوع7 20 ۲0۲ 
5 و ( 

عصلود ۷۱۵۱0 )2ع۱۱ 0ص فاصمانط ارمو ممعب۱عه وتطمصمتاملع: مطا مم صمنا2ع۹01ع1۳۷6 م۸ .2013 ریکانه ,لفط طلم5و۵ ۴6228 
ادا مصح وممممتهه لمتبهانم‌تنوه ۵۶ بالوهب نا مدوتهن رولوعط]. .م9 .رق01)) فصعماوری صمتتمصممنصا مومع 
(وتقحصحصیاک حاوتاعص طز صهتونعظ ص) .۵ 71 .صقع1 رقعع‌تتاموم۴ 

۵ 0۴ 90۵0۲7 .2008 شیک متطفاهعهام۴ مه ریگ بتطعنطع۲۱۸ 4دط۸ 24۲ ریم رتافقصطه 1 طقطی ماگ ب2زمصهمتمصهو 
۳۵۸ 8]) .12)45(:11-19 8016866 9011 20 میج ۷۷ 0۶ امتا0ل :0242 ]۱۷۲۵۲ عصتمته 40فططفع۷ مد و61 ۵2۲ ۵ صمتا0ع۲۵]1 
9 0 (2۷صصصصت ٩‏ طفتاععظ )۷1 

0 ۹۵۵ ۷۱۷۳۸ .عانااتافصا فه۸و .ع4نباع وععتعع۵ع0 9.4 و8۸ معوظ ,2015 .متام ففو 

مامصع مد ۵061 ققعهمعم ۵ ال ۷۱۵۱0 وم ۵۶ باتعصمع0عام نادمه عصتاهتمتعاظ1 ,2011 با رمتاععصلر 40ص2 بط رنه 
1 1,006هعع.[/16 ۳)6۵5://001.0۲8/10.10 .2530-2541۰ :(222)14 ۱۷۲۵0811 لهعنعم۵0۱ع۴ .عصلقومو 

۶ ۷۱۵14 240 تفه طاسمو روماطهتنه۷ ل2تامع۵ .2005 ,نک رتلل‌تهمصصها 200 و.۷.؟ مقطمظ بل ,11065 
1 1590/8010 .10 /105://001.078 . 199-207 :(62)03 012ععه۸ 12اصع1ه5 .عص2ه۲هعاه 

۶ ول ,ماجد۱ ۲۲۳۷۲۲ ]129052 عصلوتا ومتتقصصتاوع ۷۱۵۱0 م6 .2015 نا ,عطوم-60180ظ 220 ر.ل ,1۳۵۲۲۵6۵01 ریش باه17 
6 .7/1 ۲0۹://001.012/10.111) .1-16 :9 عصتوجعو م)مصصفک 0ع۵11مخ 

0 ممتاتهممم لصو ۵۶ وتانطمتنه۷ لمتاومه ملهموممطط‌اهن ,2011 .نبا بکل۷۱۵ 220 ربا ,قصفصه 1 ر.ظیت) مطعهصناطهلوم 1[ 
124-۰ :117 موم مععل۲11] 9۵11 بهلممتطاظ ‏ صصمطنمص ها اصمصمممصممه. للمو . ملمعموعالو صم.. فصمت)ه1۳0۵116 
5 105://001.0182/10.1016/[.1-01 

10۳ ۵۲۱00 0عصاصاصمی ه مه عصتونع! رومرمصه للم ۵۶ صمفتنهم‌ندمع ۸ .2000 بب۷ رمق 4ص2 ر.1 ب2عنص1 وه باه‌هالا 
199-۰ :(96)3 060۵067882 .صتقاظ ۲۴۵۸۷۵۵۵-۷۸)2822 عمط طا عامفلهع۲ ۵۶ وتمدری ۱610 لمح باتعصعل لاه ممتام0ع۳م 
۱۲۵5://001.0۲82/10.1016/50016-7061)99(00055-5 

-امصتاصهگ تما جمتاقحصطتاوع ۷1۵۱0 1166 ۳۱۵۱0-8۵216 .2019 ,.0) بلهبا 200 ,.ت) بتا2 و.م) مطلاف وک وال و.ل ,9020 رم م2 و.ل ,۷۷۵0۵ 
:(11)19 عط51ع۵ م)محصفگ بصن رممص0۳۵ باوع‌صم1 ۵۶ ممنوع؟ مصئلهه لداقدمع ما 026 (کتمق) ت۴4 سوم »م5 1۸ 
4 ۵ .2274 

1 ,0۱006۹)6۲ بلانا رقع0 ی ملز ۷۷ مطول .عافتاممتهگ آقمممصممتطظ 10۲ ومتافتاماد0ع06 .2001 .۷ وه۷نا0 20 ویک رکعا۵۵ ۷۷ 
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اعلایی و همکاران. مقایسه روش‌های زمین آمار و سنجش از دور در بر آورد ... ۵٩۹٩‏ 


11 24 مععنصا مملاداعوع۷ عصلمه ۵ وانازمانمجمع1 16 .2019 .۷۲,۲۲ رتل2اهطلوه]۱۷ 200 و.ظ رتتتحه ر.ط م221۲ مب بتاامطمه۷ 
:۱ من اماممصصمص... ۵ . لقصتام. صفتاوظ۳ .مد .. مم ...ملعم ها مها 
68 1 009://001.012/10 

۱-11610ظ)۱ جمع۱۳ع0 ومتطامممت)ه۱ع۲ لمتاهمه مه لهممجصع1 2005۰ ومسا.۷ رعصانط ۷۷ 24 ربا مطتاولا و.گ.ظ ,22100-16202 
641-۰ :(97)3 لقصدم رجمجم‌تید .ععمصا مصتفصهه مامصعد آهتامهموتهمط لهتاهمو-طونط مضه صمتامی ها واتاامامتنه۲۷ 
7 ۳۱۱0۹://01.0۲2/10.1/221081 

۴ مصتاهع1۵01 10 ومعتصاً ممتلقاهع۷6 02800۷ طمصصط0ع 0 ۸۵۵۲۷515 .2008 ,۷۷۰ ,مهن 200 رب پلالنا و. ۷ بطق ویک و۷20 و.۷ بلا2 
فصح صملاجبعوها) طاتقط . معنامورض ۵۶ . تمد تقصمتاحصعا1.. عمزر ...صرح عقمطاه . صا.. فصمتاهانمصنهعج .. ممع10۲0ظ 
6/2206 ۲۲0۵۹://01.0۲2/10.101 .1-10 :(1001 طمتامصم]ص1معی 

عصلوصعو م)مصصم؟ :ص2۵ زانط صا متبانه‌ هه 10۲ فتورولممه باتلزطاه) ناه م1 ۵11)612- ۷۲۵۱۵ :2015 رو ۷ یهععطظ 20 .یک رتقعاما70 
300-10۰ :118 سانه‌تتعض۸ 18 ۱۱ 220 کعامرصوم 00 015 220 
005://001.012/10.1016/[.001:028.2015:.6 


